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SUMMARY 


The results of an investigation of the effects of elastic deformations on the 
aerodynamic characteristics for a high-aspect-ratio, force/pressure, supercritical- 
wing model during wind-tunnel tests are presented. 

In this study a finite-element model of the wing was developed and, for condi- 
tions corresponding to wind-tunnel test points, experimental aerodynamic loads and 
theoretical aerodynamic loads were applied to the finite-element model. Comparisons 
were made between results caused by experimental data loads and analytical data 
loads, for changes in structural deflections and for changes in aerodynamic charac- 
teristics, both caused by increments in wing angle of attack and control-surface 
deflection. 

The results from this study showed good correlation between the analytical and 
experimental data. They also showed that the deformations were quite small and that 
the e>q)erimental pressure data were not significantly affected by model deformation. 


INTRODUCTION 

A series of wind-tunnel tests were conducted by the National Aeronautics and 
Space Administration (NASA) in the Langley Transonic Dynamics Tunnel to investigate 
the effects of oscillating control surfaces on unsteady aerodynamics. This series of 
tests was conducted using a model of the high-aspect-ratio (10.76), force/pressure, 
supercritical-airfoil, semispan wing that has a planform representative of current 
energy-efficient designs. The results from two of these tests are reported in refer- 
ences 1 and 2. There is concern that the model may experience elastic deformations 
that may affect the measured wind-tunnel test results, although the wing model was 
designed to be very rigid. These concerns are addressed in this paper by developing 
a finite-element structural model and combining it with an aerodynamic model to ana- 
lytically examine the model deflections caused by steady subsonic and transonic aero- 
dynamic loadings on the wing model. The purpose of the present investigation was to 
determine the deflection characteristics of the "rigid" wing model and to determine 
what effect, if any, these deflections had on the aerodynamic characteristics. 


SYMBOLS 


c . 
ra, le 

S 

C 

P 

Ac ^/Aa 


section lift coefficient 

section pitching-moment coefficient about the leading edge 
pressure coefficient 

lifting-surface pressure coefficient, C - C 

^ ^ X o p f US 

increment in section lift coefficient per change in angle of attack, deg”^ 




y spanwise coordinate, m 

z vertical coordinate, positive up, m 

a angle of attack, deg 

6 control-surface deflection, deg 

p fraction of span 

9 torsional rotation (twist), deg (see figs. 7 through 9) 


A0/Aa increment in torsional rotation per change in angle of attack 

A0/A6 increment in torsional rotation per change in control- surface deflection 

0 torsional slope along 40 percent chord, rad (see fig. 5) 

(|) bending slope along 40 percent chord, rad 
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Subscripts ; 


Is lower surface 

us upper surface 

MODELS 

Three models of the high-aspect- ratio, semispan wing are used in this study. 
They are the wind-tunnel model and two analytical models, namely, the finite-element 
structural model and the aerodynamic model. 


Wind-Tunnel Model 

. The wind-tunnel model is shown mounted in the test section of the Langley 
Transonic Dynamics Tunnel in figure 1. A sketch showing some of the geometric 
details of the model is presented in figure 2. The model has an aspect ratio of 
10.76, a leading-edge sweepback angle of 28.8°, and NASA supercritical-airfoil 
sections. Only 1 of the 10 control surfaces (number 9) was considered in this 
study. The wing model was designed and fabricated to be "rigid" by adding stiff 
boron-filament inserts to the aluminum wing-box construction as shown in figure 3. 
The model was attached to a five-component force balance which in turn was attached 
to a remotely controlled turntable mounted in the wind-tunnel wall. A complete 
description of the model is given in references 1 and 2. 


Finite-Element Structural Model 

An elastic-axis, finite-element structural model was constructed to represent 
the wind-tunnel semispan wing model. The finite-element model was developed using a 
system of computer programs called SPAR (now called EAL, Engineering Analysis Lan- 
guage) . SPAR was designed primarily to perform stress, buckling, and vibrational 
analyses of linear finite-element systems (ref. 3) . Figure 4 shows a sketch of the 
finite-element model used in this study. The model of the wing planform consists of 
62 beam elements, 2 quadrilateral plates, and 1 zero-length beam spring element. 

These elements are positioned by nodal points which are listed in table 1 . Twenty of 
the beam elements are along the elastic axis, the 40-percent-chord line. The stiff- 
nesses of the 40-percent-chord beam elements were determined from laboratory mea- 
surements on the wind-tunnel model. The remaining beam elements were assigned a 
rigid stiffness (El = 5.74 GJ = 2.87 t^^)^^nd were oriented perpendicular to 

the elastic axis so that deflections at points off the elastic axis could be 
determined. The location of the nodal points for the beam elements along the elastic 
axis were selected to correspond to pressure-measurement stations and the locations 
midway between the pressure-measurement stations. The two quadrilateral plates model 
the control surface and were also assigned a rigid stiffness. The zero-length beam 
spring element represented the root- attachment flexibility, torsion and bending, of 
the model mounting system. The stiffnesses used for this element were also deter- 
mined from laboratory measurements. 

Measured and calculated structural data are presented in figure 5. These 
results were obtained by applying concentrated loads near the tip. A vertical load 
of 444.8 N (100 Ibf), gave the measured and calculated bending-slope results shown in 
figure 5(a). A torsional load of 11.298 N-m (100 Ibf-in.) gave the measured and 
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calculated torsional-slope results shown in figure 5(b) . These results, determined 
along the elastic axis, indicate that the structural modeling of the wind-tunnel 
model was very accurate. 


Aerodynamic Model 

The aerodynamic model was generated by using linear subsonic lifting-surface 
theory as implemented through the RHOIV digital computer program (refs. 4 and 5) . 
RHOIV was developed to calculate unsteady aerodynamic loadings caused by motions of 
lifting surfaces with leading-edge and trailing-edge control surfaces. RHOIV was 
selected to model the aerodynamics in this study because of its ability to incorpo- 
rate control surfaces in its calculations and because of the ease with which it can 
be used. 

For this study, steady-state aerodynamic loadings were determined from RHOIV 
when the reduced frequency was set to zero. A sketch of the RHOIV aerodynamic model 
is shown in figure 6. This figure shows the points used to define the planform 
shape. The locations of the collocation points used to calculate the pressures were 
determined through a default option available in RHOIV. This program treats the 
control surface and the wing separately and accounts for aerodynamic flow singulari- 
ties of the control-surface hinge line. The aerodynamic model is a flat plate with 
the planform shape of the wind-tunnel model including one outboard trailing-edge 
control surface. When displacement points were used to represent a deformed wing 
shape, the location points were the same as the nodal points in the finite-element 
structural model. 


ANALYSES 

The finite-element structural model was combined with the RHOIV aerodynamic 
model to study the elastic deformation effects on the aerodynamic loading of the wing 
model at M = 0.60 and M = 0.78. The flow about the model at M = 0.60 was a 
subsonic test condition at which linear aerodynamic theory would be expected to be 
accurate. This aided in the comparison of experimental and analytical results. Mach 
0.78, a transonic test condition, is the design Mach number for the wing. Analysis 
at both Mach numbers determined the deformations of the wing model due to a change in 
angle of attack and the aerodynamic characteristics due to these deformations. This 
analysis was also done at M = 0.60 for an increment of control-surface deflections. 
The aerodynamic load conditions were of two forms. One analysis used selected exper- 
imental pressure results from reference 2, and the other analysis used aerodynamic 
loads calculated by using the RHOIV aerodynamic model. 


Analysis With Experimental Loads 

The model deformations calculated from the experimentally determined wing- 
surface loadings were obtained by the following procedure. For the conditions of 
M = 0.60 and q = 3.033 kPa, the pressure distributions from reference 2 (points 159 
and 214) were used. One set of data was for a = 2.85°, and the other was for 
a= 0.0°. The data were numerically integrated chordwise to obtain concentrated 
forces and moments at appropriate points along the 40-percent-chord line (elastic 
axis) of the wind-tunnel model. The concentrated loads were representative of a 1° 
change in angle of attack. The concentrated forces and moments were then applied to 
the corresponding nodal points of the structural model, and the resulting vertical 
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displacements and rotations along the elastic axis were calculated. This procedure 
was also for M = 0.78 (q = 3.84 kPa, a = 2.05° (point 63), and a = 0.0° 

(point 55)) and for M = 0.60 (q = 3.052 kPa, a= 2.85°, 6 = -6° (point 184), 
and 6= -2° (point 178)). (See ref. 2.) 


Analysis With Theoretical Loads 

The model deformations calculated from theoretically determined wing-surface 
loadings were obtained by using an iterative procedure. The RHOIV aerodynamic pro- 
gram was used to determine wing-surface pressure distributions for the same Mach 
numbers as the experimental loads analysis. However, at M = 0.60 and M = 0.78, 
the wing angle of attack was initially defined at 1°. Also, at M = 0.60 the wing 
angle of attack was initially defined at 0°, and the control surface was set at 1° 
deflection. The calculated pressure distributions were then converted to concen- 
trated forces and moments along the elastic axis and applied to the finite-element 
model. The resulting deformation formed a new wing shape. The new wing shape was 
then used with the aerodynamic model to determine a new aerodynamic loading, which 
was applied to the finite-element model. The process was repeated until there was 
essentially no difference in the pressure distributions for successive iterations. 
In general, the process converged rapidly, requiring only two iterations. 


DISCUSSION OF RESULTS 

This investigation was conducted to determine the deflection characteristics of 
the wing model and to determine if these deflections had any effect on the aerody- 
namic characteristics. The analysis was conducted by using a finite-element struc- 
tural model and an aerodynamic model, which were discussed previously. The results 
of the analysis are presented in the following sections. 


Structural Results 

A method of determining the rigidity of the wing model is to examine the defor- 
mations caused by aerodynamic loadings on the wing. As discussed in the section 
entitled "Analyses," two types of loads, experimental and theoretical, were applied 
to an accurate structural finite-element model of the wind-tunnel wing model. The 
deformations resulting from these loadings are shown in figures 7 through 9. Fig- 
ure 7 shows the bending and torsional deflections for M = 0.60 at an incremental 
angle of attack. Similar results are shown in figure 8 for M = 0.78. In figure 9, 
the bending and torsional deflections are shown for M = 0.60 at an incremental 
control-surface deflection. The deformations shown in figures 7, 8, and 9(a) are, in 
general, similar in shape and value, although the theoretical and experimental load 
results are not in exact agreement. However, figure 9(b) shows a large discrepancy 
between the two load cases, the cause of which has not yet been determined. In all 
the figures, the differences between the theoretically and experimentally determined 
deformations are numerically small, but more importantly the overall deformations are 
considered very small. For example, at M = 0.78 for an incremental angle of attack 
the change in bending deflection is determined to be no more than 1.8 mm per angle of 
attack, and the change in rotation is not more than 0.025° per angle of attack at the 
tip where the greatest deflections occur. (See figs. 8(a) and 8(b).) The shape of 
the deflections along the wing span are similar to previous studies of wing deforma- 
tions. (See ref. 6.) An attempt was made to measure actual deformations as they 
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occurred during the wind-tunnel tests by using a cathetometer . Unfortunately, the 
resolution of the cathetometer was not high enough to give accurate results. 


Aerodynamic Results 

Another method used to determine the rigidity of the wing model is to examine 
the effects the deformations had on some of the aerodynamic characteristics of the 
wing model. The results of this study are shown in figures 10 through 15. 

As discussed in the section entitled "Analyses," an aerodynamic model was used 
with the RHOIV program to determine theoretically the pressure distribution over the 
wing model. Figures 10 through 12 show the pressure distributions at one of the 
semispan locations (ti = 0.71) for the different Mach numbers with the incremental 
angle of attack and control-surface deflections for the experimental and theoretical 
loads. At M = 0.60 for an incremental angle of attack (fig. 10), three sets of 
data are shown, the experimental data, the data from the aerodynamic model without 
deformations, and the data from the aerodynamic model with the final deformations of 
the iterative procedure. As can be seen, there is fairly close agreement between the 
experimental data and the two sets of theoretical data. The close agreement between 
the two sets of theoretical data shows that the deformations created by the aerody- 
namic loads do not significantly affect the pressures at this semispan location of 
the wing. Similar results at an incremental angle of attack for M = 0.78 and at an 
incremental control-surface deflection for M = 0.60 are shown in figures 11 
and 12, respectively. In all three figures there is some disagreement between the 
experimentally and theoretically determined data. Some of the discrepancy could be 
due to slight scatter, but another source of slight error comes from the RHOIV aero- 
dynamic program. As the program was initially used, the dynamic pressure distribu- 
tions at transonic Mach numbers and over control-surface hinge lines did not give as 
accurate results as are now available (ref. 5); however, the results achieved were 
considered adequate for this study. 

Figures 13 through 15 show the change in section lift coefficient and section 
pitching-moment coefficient over the wing span at M = 0.60 for an incremental angle 
of attack, M = 0.78 for an incremental angle of attack, and M = 0.60 for an 
incremental control-surface deflection, respectively. The experimental and theoreti- 
cal pressure data were integrated over the wing chords to give these results, which 
were used with the finite-element structural model to determine the deformed shape of 
the wing model. At M = 0.60 for an incremental angle of attack (fig. 13), there is 
fairly close agreement for the section lift coefficient over the wing span. There is 
much more scatter in figure 13(b) for the section pitching-moment coefficient over 
the wing span. Some of this scatter, as in figures 14 and 15, may have resulted from 
some nonlinearity between angles of attack and control-surface deflections in deter- 
mining the lifts and moments for incremental values. Figure 14 shows the section 
lift and pitching-moment coefficients at M = 0.78 for an incremental angle of 
attack. In figures 14(a) and 14(b), there is still some disagreement between the 
experimental and theoretical cases, although there is better agreement for the moment 
at M = 0.78 than at M = 0.60. Also, the two sets of theoretical data compare well 
for the section lift and pitching-moment coefficients. Figure 15 shows the section 
lift and pitching-moment coefficients at M = 0.60 for an incremental control- 
surface deflection. These results are similar to those already presented. However, 
there are only two sets of data in figure 15(b). The theoretical data with the 
deformation effects for an incremental section pitching-moment coefficient were not 
available. 
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SUMMARY OF RESULTS 


The concerns that the wind-tunnel model would experience deformations and that 
these deformations would affect the aerodynamic characteristics of the wing have been 
addressed and examined through two methods of analysis. The results of this study 
indicate the following: 

1. The deformations of the wing model are very small. That is, the maximum 

deflections are at M = 0.78, where the increment in bending deflection per 
change in angle of attack is 1.61 mm/deg and the increment in torsional 
rotation per change in angle of attack is SOBS, 0 'OZ.O^S, 

2 . There is good agreement between the theoretical and experimental pressure 

data. 

3. The iterative procedure used to calculate the theoretical aerodynamic charac- 

teristics shows that there is little difference between the results for the 
deflected and nondeflected wing models. 

4. The good agreement between the theoretical and experimental aerodynamic data 

and the close agreement between the theoretical aerodynamic results for the 
deflected and nondeflected model shapes shows that the deformations do not 
significantly affect the aerodynamic characteristics of the wing model. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
April 14, 1982 
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Figure 6.- Sketch of RHOIV aerodynamic model. 
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Figure 10.- Incremental lifting-surface pressure distribution at n = 
and M = 0.60 for an incremental angle of attack. 
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Figure 14.- Spanwise aerodynamic characteristics of wing model at 
M = 0.78 for an incremental angle of attack. 
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(a) Incremental section lift coefficient. 



(b) Incremental section pitching-moment coefficient. 


Figure 15.- Spanwise aerodynamic characteristics of wing model at M 
for an incremental control-surface deflection. 
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